The Kondo Effect in Ferromagnetic Atomic Contacts 
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Iron, cobalt and nickel are archetypal ferromagnetic metals. In bulk, electronic conduction in 
these materials takes place mainly through the s and p electrons, whereas the magnetic moments 
are mostly in the narrow d-electron bands, where they tend to align. This general picture may 
change at the nanoscale because electrons at the surfaces of materials experience interactions that 
differ from those in the bulk. Here we show direct evidence for such changes: electronic transport 
in atomic-scale contacts of pure ferromagnets (iron, cobalt and nickel), despite their strong bulk 
ferromagnetism, unexpectedly reveal Kondo physics, that is, the screening of local magnetic mo- 
ments by the conduction electrons below a characteristic temperature^. The Kondo effect creates 
a sharp resonance at the Fermi energy, affecting the electrical properties of the system;this ap- 
pears as a Fano-Kondo resonance^ in the conductance characteristics as observed in other artificial 
nanostructure o^'"^'^'^'^'^'^'"^'^'"^"^ . The study of hundreds of contacts shows material-dependent lognor- 
mal distributions of the resonance width that arise naturally from Kondo theoryi^. These resonances 
broaden and disappear with increasing temperature, also as in standard Kondo systems^ '^'^'^ . Our 
observations, supported by calculations, imply that coordination changes can significantly modify 
magnetism at the nanoscale. Therefore, in addition to standard micromagnetic physics, strong elec- 
tronic correlations along with atomic-scale geometry need to be considered when investigating the 
magnetic properties of magnetic nanostructures. 



Atomic-scale contacts can be fabricated by techniques 
such as scanning tunnelling microscopyi^ or the use of 
electromigrated break junctions (EB Js)i^, where the size 
of a macroscopic contact between two leads is reduced 
until they are in contact through only a few atoms and, 
eventually, through only one. The conductance of metal- 
lic monatomic contacts is known to be around 2Go, where 
Go = e^//i is the spin-resolved quantum of conductance^'^ 
(e being the elementary charge and h Planck's constant). 
To identify the atomic contacts, histograms are con- 
structed from the evolution of the conductance recorded 
during the breaking of different contacts (Fig. la, b). 
The position of the first peak of these histograms is iden- 
tified as the conductance of the monatomic contact. For 
iron, cobalt and nickel, the conductance is larger than 
2Go owing to the contribution of the sp and d orbitals to 
the transmissio n^^i^^i^^ . 

We have studied the low-temperature conductance 
characteristics of hundreds of atomic-scale contacts of 
the three transition-metal ferromagnets iron, cobalt and 
nickel using a home-built STM. More than the 80% 
of the differential conductance {dl/dV) curves at the 
monatomic contact show peaks or dips around zero bias 
such as those shown in Fig. Ic, which are very similar to 
those observed in STM spectroscopy of single magnetic 
adatoms on non-magnetic surface o^'^'^i^^ . Thus, as in the 
case of these Kondo systems, we can also fit our dl/dV 
curves to the sum of a flat component, go, and a Fano-like 
resonance that typically amounts for 10% of the signal: 
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Here e = {eV — €s)/kBTK is the bias shifted with re- 
spect to the centre of the resonance, eg, and normalized 
by the natural width of the resonance, ksTn', Tk is the 
Kondo temperature; q is the dimensionless Fano param- 
eter that determines de simmetry of the curve; A is the 
amplitude of the feature; and ks is Boltzmann's con- 
stant. 

It is clear from the dl/dV characteristics in Fig. Ic 
that the width of the Fano feature is different for each 
material. This is confirmed by statistical analysis of the 
data. Figure 2 shows histograms of Tk obtained from 
fitting our conductance curves for hundreds of different 
iron, cobalt and nickel monatomic contacts. Notably, 
these histograms follow a log-normal distribution; that 
is, the natural logarithm of Tk is normally distributed. 
Because many different atomic configurations result in 
monatomic contacts, their electronic properties, such as 
conductance (Fig. lb), density of states and the asso- 
ciated energy scales, are expected to be normally dis- 
tributed. Instead, a normal distribution of log(TK) can 
only be understood if Tk can be expressed as the expo- 
nential of normally distributed quantities. The Kondo 
model^^ naturally relates Tk to the exponential of the 
typical energy scales in the problem. 

Fitting the histograms to a log-normal distribution 
yields most frequent values for the resonance widths in 
the different materials corresponding to Tk—90 K, 120 K 
and 280 K for iron, cobalt and nickel, respectively, follow- 
ing the same trend (Tj^" < T^° < T^') for these chemi- 
cal species when deposited as adatoms on non-magnetic 
surfaces^. In simple terms, the Kondo temperature de- 
creases as the size of the screened magnetic moment in- 
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FIG. 1: Conductance of a monatomic contact, a, Example of a trace where we record the conductance while stretching a nickel 
wire using a scanning tunnelling microscope (STM) at 4.2 K. Inset, model of a monatomic contact, b, Conductance histograms 
constructed for iron, cobalt and nickel from thousands of such traces. The position of the first peak of in each histogram 
corresponds to the conductance of the monatomic contact, c. Differential conductance curves recorded at the monatomic 
contact as a function of the applied voltage. A characteristic resonance appears at small bias that fits the Fano line shape. All 
possible symmetries are found in the spectroscopy of iron, cobalt and nickel contacts, and the width of the resonance is the 
main difference between the spectra of the three materials. This width is proportional to the Kondo temperature. 
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FIG. 2: Histograms of inferred Kondo temperatures for iron, 
cobalt and nickel. The histograms are constructed from more 
than 200 fittings and normalized to the total number of curves 
fitted. The continuous lines show the fits of the data to log- 
normal distributions of Tk with a different most probable 
value for each material. 

creases, as we change from nickel to cobalt to iron. The 
Kondo temperature of cobalt nanocontacts is very sim- 
ilar to the one reported in reiM- for cobalt in Cu(lOO) 
probed using a STM in the contact regime. 

In addition to the statistical analysis described above, 
we measured the temperature evolution of the Kondo fea- 
tures for a single contact. To preserve the atomic stabil- 
ity of the junction while changing the temperature^^, we 
used an EBJ such as the one shown in the inset of Fig. 
3b. It is important to note that in this alternative imple- 
mentation we observe exactly the same Fano resonances, 
with the same distribution of Kondo temperatures, as we 



obtained in the STM experiments using nickel and cobalt 
samples. In all the cases, we see a reduction of the ampli- 
tude of the Fano features as the temperature increases, 
as shown in Fig. 3a for the case of cobalt. In Fig. 3b 
we plot A{T), as defined in equation (1), as a function 
of T on a logarithmic scale. The curve so obtained has a 
low-temperature plateau followed by a linear decay, very 
similar to that of quantum dots and molecules in the 
Kondo regime^'^'^S. 

In summary, our atomic contacts fabricated with two 
different methods show the same dl/dV curves (Fig. Ic), 
the same chemical trends (Fig. 2) and the same tempera- 
ture evolution (Fig. 3a, b) as other standard, chemically 
inhomogeneous Kondo systems. This indicates that the 
contact atom(s) in nanocontacts of iron, cobalt and nickel 
are in the Kondo regime. This is unexpected for two 
reasons. First, the Kondo effect has always been associ- 
ated with chemically inhomogeneous systems containing 
at least two kinds of atom: those where the localized level 
resides and those providing the itinerant electrons. Here 
the same chemical species hosts both the itinerant states 
and the local magnetic moments. Second, in most re- 
ports on the Kondo effect, the itinerant electrons are not 
spin-polarized, as a large spin polarization is expected 
to destroy the effect, or a carbon nanotube^i contacted 
with ferromagnetic leads has been justified by the op- 
positely directed spin polarizations in the electrodes, in 
agreement with refsj^i^i^i. Although the situation re- 
garding the magnetization of the leads might be similar 
in our system, we argue that the Kondo effect is still 
possible even if there is no domain wall pinned in the 
contact. 

Being in the Kondo regime implies that the atoms at 
the contact must host, at least, a localized d-electron level 
whose magnetic moment is screened as a result of antifer- 
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FIG. 3: Evolution of the Fano resonances with increasing tem- 
perature, a, Characteristics of differential conductance ver- 
sus bias voltage for a cobalt atomic contact, showing how the 
Fano resonance disappears as the temperature is increased, 
b, The amplitude of the Kondo resonance from a decreases 
logarithmically with the temperature. Inset, an artificially 
coloured example of a lithographic device for the EBJ exper- 
iments (before electromigration) ; the cobalt junction is blue. 
The atomic contacts created by this method are suitable for 
studying the temperature dependence of the Fano resonance. 



romagnetic coupling to the sp conduction electrons. We 
can show how local moment formation and antiferromag- 
netic coupling occur in nickel nanocontacts in the follow- 
ing way (the cases of iron and cobalt can be understood 
on similar grounds) . The mean- field solution of the An- 
derson model^^, which describes a localized d level, with 
energy e^, on-site repulsion U and hybridization with the 
itinerant sp electrons, T^p-d, defines the conditions for 
the formation of a local moment in such a d level. The 
model can also be used to derive the antiferromagnetic 
sd exchange coupling, jff^ (ref.^^), which arises from 
the sp — d hybridization term. Here we used the local 
spin-density approximation (LSDA) to density functional 
theory and its generalization LSDA-I-U as the mean field 
to determine the formation of a local moment and its 
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FIG. 4: Electronic structure for a nickel chain and a nickel 
nanocontact. a, LSDA minority-spin energy bands for a 
monostrand nickel chain, where the wavevector k runs over 
the first Brillouin zone. The bands are labelled according to 
their symmetry group: the El and E2 bands (d-band like, 
doubly degenerate and decoupled from the s band) and the 
Al or d^^2_^2 band (hybridized with the a band). The lat- 
tice constant is 2.09 A . b. As in a, but for majority-spin 
electrons. The yellow background indicates occupied states, 
c, LSDA-I-U minority-spin results for the total DOS and the 
DOS projected on the d^2_y2 orbital for a tip atom of a nickel 
nanocontact, with (7= 3 eV. d, As in c, but for majority-spin 
electrons. 



antiferromagnetic coupling to the sp carriers. We con- 
sider both monostrand chains and nanocontacts. Their 
smaller atomic coordination, compared with that of the 
bulk, results in a stronger electronic localization and a 
larger magnetic moment per atom (values of 1.17 Bohr 
magnetons, compared with 0.6 Bohr magnetons in the 
bulk, were obtained using LSDA). This favours the ap- 
pearance of the Kondo effect. In Fig. 4a, b, we show the 
LSDA bands obtained for the nickel chain. Out of the 
six minority-spin bands crossing the Fermi level, the two 
degenerate E2 bands are the narrowest, hosting highly 
localized electrons. These are the bands that are less 
well described by the LSDA because of the inherent self- 
interaction problem of this approximation. 

We modelled actual nanocontacts with two identical 
pyramids facing each other in the [001] direction. As pre- 
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viously done for chains^^, and to avoid the self- interaction 
problem, we computed the electronic structure of the 
nanocontact using LSDA+U. In Fig. 4c, d, we show the 
total density of states (DOS) projected on a tip atom 
and the DOS projected on the corresponding d^2_y2 oi- 
bitals for both the minority (Fig. 4c) and the major- 
ity (Fig. 4d) spins (here we have set U=3 eV). In gen- 
eral, solutions obtained for values of U ranging from 3 
to 5 eV also show that the dx2-y2 orbital, forming one 
of the E2 bands in chains, hosts an integer local mag- 
netic moment. Different geometries may favour the for- 
mation of the local moment in other strongly localized 
orbitals. The hybridization of the dj.2_y2 orbital with 
the surrounding sp orbitals results in antiferromagnetic 
kinetic exchange^^. We estimate that J^^^ — 1 eV, 
taking ~ 5 eV, |I4p-dP ~ 2 eV and Ueff ^ S eV 
from our LSDA-I-U calculations, where Ueff is the spin 
splitting of the dx2 -y2 orbital (Supplementary Informa- 
tion). We note that, in contrast to the contact atom(s), 
the intra-atomic sp — d hybridization vanishes for bulk 
atoms owing to the very small anisotropy of the crystal 
environment. This makes the antiferromagnetic coupling 
J^^d larger in the contact than in the bulk. This cou- 
pling competes with the ferromagnetic coupling JfJ^I^^ 
which we estimate from the splitting of the sp band at 
the boundary of the Brillouin zone for chains to be ~0.2 
eV (Fig. 4a, b). Thus, an overall antiferromagnetic cou- 
pling between d electrons and sp conduction electrons is 
possible in nickel nanocontacts. Additionally, the local 
moment, m, responsible for the Kondo effect is subject 
to the ferromagnetic coupling Jdd to the neighbouring 
atoms, and this interaction also competes with the anti- 
ferromagnetic sp — d coupling: 
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Here the rrii are the local moments of the neighbouring 
atoms and Sg is the spin of the s electrons. In ref^SS,^ the 
coupling Jdd was calculated for iron, cobalt and nickel by 
implementing the magnetic force theorem with a LSDA 
ground state. This method yields values for the spin- 
wave dispersion of the materials that compare well with 
experiment^, and gives Jdd=19 meV, 15 meV and 2.7 
meV for iron, cobalt and nickel, respectively. Thus, Jdd 
is significantly smaller than Jfp^d- 

The Kondo effect in nanocontacts is favoured by three 
factors. First, a local moment forms in the contact 
atoms because of their smaller coordination. Second, 
the reduced symmetry of the contact, compared with 
the bulk, enhances the intra-atomic contribution to the 
sp—d hybridization and, thus, the antiferromagnetic cou- 
pling JfJLd ■ Third, the smaller coordination also re- 
duces the influence of the direct ferromagnetic dd cou- 
pling with neighbouring atoms. As a result, the local 
moment formed in the contact is antiferromagnetically 
coupled to the sp itinerant electrons and results in the 



Kondo effect in this system, in contradiction to conven- 
tional wisdom. 
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APPENDIX A: METHODS 

For the statistical results on the Kondo parameters, 
we used a home-made STM at 4.2K to fabricate the con- 
tacts by indentation between two pieces of metal^^. For 
the temperature-dependence measurements, the contacts 
were produced by the controlled electromigration at 4.2K 
of 100-nm-wide junctions fabricated by electron beam 
lithography^^. In both cases, the spectroscopic curves 
were obtained by the addition of an a.c. voltage with a 
peak-power amplitude of ImV and a frequency of 1 kHz 
to the d.c. bias voltage to allow the lock-in detection of 
the differential conductance. 

a. Fabrication of atomic-scale contacts by scanning 
tunnelling microscopy. Two pieces of metal wire of 0.1- 
mm diameter were cleaned and sonicated in acetone and 
isopropanol. These pieces were mounted in a home-built 
STM. The set-up was pumped down to high vacuum and 
immersed in a liquid helium bath until the temperature 
reached 4.2 K. The two pieces of wire were brought into 
contact and then pulled apart until the contact was of 
atomic dimension and, finally, until only one atom formed 
the contact^^. Strong indentation before the formation of 
every contact ensured the cleanliness of the atomic con- 
tacts. The histograms obtained by this technique (Fig. 
la) show similar results to the ones obtained by the me- 
chanically controlled break junction technique, where the 
surfaces brought into contact are created under cryogenic 
conditions. 

b. Fabrication of atomic-scale contacts using EBJs. 
As a first step, a small junction about 100 nm wide 
was fabricated from cobalt using electron beam lithog- 
raphy and electron beam evaporation over a silicon diox- 
ide substrate. To make the junction suitable for con- 
tact by macroscopic probes, two gold electrodes are de- 
posited over the edges of the junction in a second lithog- 
raphy step, following the procedure described in ref.^». 
The samples are then placed in a probe station that 
was pumped down and immersed in a liquid helium bath 
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until the sample reached a temperature close to 4.2 K. 
Under these conditions, the controlled electromigration 
proces o^^'^^ was performed, decreasing the size of the 
jimction to the atomic scale. 



APPENDIX B: SUPPLEMENTARY 
INFORMATION ON LSDA AND LSDA+U 
RESULTS 

The electronic structure of a Ni nanocontact model 
as that shown in Fig. 5 has been calculated with den- 
sity functional theory in the local spin-density approx- 
imation (LSDA)"^°using our established ab-initio trans- 
port methodologyi^. Fig. 5a shows the DOS and par- 
tial d DOS for one of the tip atoms of the nanocontact. 
The majority-spin d-levels are completely filled while the 
minority-spin d-levels are partially occupied with an av- 
erage occupation of 0.8. This results in a net magnetic 
moment of ~1 for the tip atom similar to the magnetic 
moment of the one-dimensional chain. 

It is, however, well known that LSDA suffers from the 
so-called self-interaction problem which blue-shifts occu- 
pied levels. Here the spurious self-interaction brings all 
the minority d-orbitals up to the Fermi level leading to a 
partial and almost equal occupation of 0.8 of all five d- 
orbitals. The self-interaction error of LSDA is corrected, 
e.g., in the LSDA-hU method^!. In the LSDA-^U method 
the effective LSDA Kohn-Sham potential of the strongly 
interacting d-electrons is corrected by adding a Hartree- 
Fock term for an on-site Coulomb repulsion U and ex- 
change interaction J: 

{ia\VLSDA+u\ic^) ^ UiNd-n'^)-J{NS~n'^)-Edc (Bl) 

Here, Nd is the total occupation of all the d-levels of an 
atom while nf is the occupation of an individual d-level 
per spin a, and is the total occupation of d-levels per 
spin a. Edc accounts for the fact that the Coulomb repul- 
sion and exchange interaction have already been taken 
into account in some way at the LSDA level, and thus 
have to be subtracted from the effective LSDA Kohn- 
Sham potential in order to avoid double-counting^. 

Due to the low coordination of the two tip atoms 
(compared to the rest of the nanocontact and to bulk 
atoms) Coulomb interaction effects are expected to be 
much stronger for the tip atoms than for the rest of the 
nanocontact. Hence, we only treat the two tip atoms of 



the nanocontact at the LSDA-I-U level while the rest of 
the nanocontact is treated at the LSDA level. For bulk Ni 
a widely accepted value is U =3 eV for the direct Coulomb 
repulsion and J=leV for the exchange interaction-^. Due 
to the lower coordination of the tip atoms with respect 
to bulk atoms, the screening of the Coulomb repulsion 
should be lower than in bulk. Therefore the Coulomb 
repulsion is expected to be bigger than [7=3 eV for the 
tip atoms. 
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\tsd\ 


\tpd 1 


Ed 


ea + U 


jAF 

^ sp — d 


3 eV 


0.3 eV 


0.6 eV 


-4.5 eV 


1.5 eV 


1.6eV 


5 eV 


0.3 eV 


0.6 eV 


-5 eV 


3 eV 


l.OeV 



TABLE I: Table of hoppings and excitation energies extracted 
from an LSDA-f U calculation with U— 3 eV and U=5 eV, and 
J=l eV for all d-orbitals on the two tip atoms, as well as the 
resulting antiferromagnetic exchange coupling jfJLd- 



Indeed, as can be seen from Figs. 5a, b, for U >3 eV 
four of the minority d-orbitals are now well below the 
Fermi energy and completely occupied, while the minor- 
ity d^2_j,2 orbital is completely empty. Thus the magnetic 
moment of the tip atom is now entirely carried by the 
c?2.2_j^2-level. This shows that a local moment can form 
on the tip atom which is carried by a single d-orbital. 
The screening of the magnetic moment of this c?-orbital 
by the sp-conduction electrons can thus give rise to a 
Kondo resonance at the Fermi level. 

Now we estimate the exchange coupling J^^^ of the 
1^3.2 _j,2 orbital of the tip atoms from the hoppings and 
excitation energies of the effective Hamiltonian of the 
LSDA-f U calculation. The da;2_j,2-orbital of the tip atom 
of one of the pyramids is coupled via hoppings tsd and 
tpd to the s and p^-orbitals of the four atoms next to the 
tip atom in the same pyramid as detailed in Tab. 1. The 
total antiferromagnetic exchange coupling is given by: 

where Vsp-d is the effective total hybridization with 
the sp-channels: 

=4|t,dP+4|ipdp (B3) 
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